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I .  INTRODUCTION 


The  properties  of  the  propellant  gas  near  the  muzzle  during  and  shortly 
after  shot  ejection  are  of  interest  since  they  determine  the  magnitude  of  the 
muzzle  blast  overpressure.  One  important  property  is  the  temperature.  Schmidt, 
Gion,  and  Shear  designed  an  apparatus  to  measure  the  sound  speed  and  thereby 
calculate  a  temperature  by  assuming  an  equation-of-state  for  the  propellant 
gas.  Schmidt,  et  al 1  measured  the  propellant  gas  pressures  as  a  function  of 
time  at  different  locations  near  the  muzzle.  For  slower  moving  projectiles, 
the  propellant  gas  behind  the  projectile  has  subsonic  flow  conditions  immedi¬ 
ately  before  shot  ejection.  After  the  projectile  uncorks  from  the  muzzle,  an 
expansion  wave  propagates  into  the  gun  tube.  They  approximate  this  expansion 
wave  as  being  one-dimensional  and  thereby  deduce  the  sound  speed.  Further 
approximations  may  be  made  to  calculate  the  propellant  temperature. 


The  complexity  of  the  starting  expansion  wave  might  contribute  toward  the 
disagreement  in  temperature  between  experimental  results  and  predictions  from 
interior  ballistics  theory."  Immediately  after  shot  ejection,  the  expan¬ 
sion  wave  entering  the  gun  tube  is  two-dimensional.  Gough2 3 4 5  contends  that  the 
initial  two-dimensional  nature  of  this  wave  might  mask  any  effects  produced 
by  the  initial  chemically  frozen  wave  propagating  into  the  tube.  He  also 
concludes  that  The  geometrical  dispersion  due  to  the  two-dimensional  geome¬ 
try  of  the  unloading  wave  is  too  complicated  to  quantify  as  a  source  of 
error  without  further  study."  With  the  NOVA  simulation  code,  Gough  also 
found  that  the  temperature  of  the  propellant  gas  near  the  base  of  the  pro¬ 
jectile  was  markedly  higher.  Additionally,  Celmins B  indicates  that  the  usual 
assumption  of  3  very  slow  variation  in  propellant  gas  temperature  with  dis¬ 
tance  from  the  muzzle  might  not  be  valid. 


1.  a.  M.  ochmidt,  E.  J.  Gion,  and  D.  D.  Shear,  "Acoustic  Thermometric 
Measurements  of  Propellant  Gas  Temperatures  in  Guns, "  U.  S.  Army  Ballistic 
Research  Laboratory ,  Aberdeen  Proving  Ground,  MB,  BRL  Report  BRL-R-1919 
August  1976  (ad  A030359) 

2.  P.  S.  Gough,  "On  the  Determination  of  the  Gas  Temperature  from  the  Velocity 
of  the  Muzzle  Rarefaction  Wave,"  U.  S.  Army  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  MB,  BRL  Contractor  Report  BRL-CR-00504 

February  1983  (AD  Al 254 7-9)  .  J 

3.  A.  Celmins,  "Theoretical  Accuracy  of  Acoustic  Gas  Temperature  Measurements 
in  Guns, "  U.  S.  Army  Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground, 
MB,  BRL  Technical  Report  ARBRL-TR-02134,  January  1979  (AD  A068461). 

4.  A.  Celmins,  "Theoretical  Basis  of  the  Recoilless  Rifle  Interior  Ballistics 
Code  RECRIF,"  U.  S.  Army  Ballistic  Research  Laboratory ,  Aberdeen  Proving 
Ground,  MB,  BRL  Report  BRL-R-1931 ,  September  1976  (AD  B013832L) . 

5.  P.  G.  Baer  and  J.  M.  Frankie,  "The  Simulation  of  Interior  Ballistic 
Performance  of  Guns  by Digital  Computer  Program,"  U.  S.  Army  Ballistic 
Research  Laboratory ,  Aberdeen  Proving  Ground,  MB,  BRL  Report  BRL-R-1183 
1962  (AD  299980). 
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-t  is  the  objective  of  this  work  to  examine  the  one-dimensional  assump- 
lon  and  the  flow  details  near  the  muzzle  after  shot  ejection.  The  time  for 
transition  from  two-dimensional  to  one-dimensional  flow  is  ascertained  and  the 
cttects  upon  t)ie  acoustic  thermometry  measurements  and  flow  exterior  to  the 
gun  barrel  muzple  is  discussed.  A  gun-blast  numerical-simulation  code6’7 
modified  for  the  particular  problem  is  used  to  calculate  the  flow  properties 
both  inside  and  outside  the  gun  barrel.  The  results  of  the  simulation  can 
be  compared  with  an  experiment  conducted  at  BRL. 8  To  examine  the  transition 
rom  a  two-dimensional  flow  to  the  one-dimensional  rarefaction  wave,  Mach 
contours  are  used  to  describe  the  different  stages  of  the  gun  exhaust  cycle. 


II.  EXPERIMENT 

In  acoustic  thermometry,  the  velocity  of  propagation  of  an  acoustic  wave 
is  found  by  measuring  the  time  for  it  to  travel  a  known  distance.  For  this 
experiment , the  disturbance  is  generated  when  the  projectile  uncorks  and  an 
expansion  wave  travels  into  the  bore  at  a  velocity  equal  to  the  propellant 
gas  speed  of  sound  minus  the  local  outflow  velocity  of  the  propellant  gases. 
An  in-bore  wave  diagram  is  shown  in  Figure  1.  Five  pressure  transducers, 
mounted  flush  in  the  tube  wall  near  the  muzzle,  record  the  passage  of  the 
expansion  wave.  A  pressure-time  curve  generated  by  the  experiment  is  shown 
in  Figure  2.  The  rise  in  the  first  part  of  the  pressure  curve  is  caused  by 
air  being  compjressed  ahead  of  the  projectile  as  it  travels  down  the  tube. 

The  rapid  rise  in  pressure  is  associated  with  the  passage  of  the  projectile 
over  the  pressure  transducer.  There  is  then  a  decay  caused  by  in-bore  waves 
reflecting  from  the  breech  and  the  base  of  the  projectile.  After  the  pro¬ 
jectile  clears  the  muzzle,  the  expansion  wave  travels  back  toward  the  trans¬ 
ducer  and,  as  it  passes,  the  recorded  pressure  declines  at  a  faster  rate.  The 
velocity  of  the  propellant  gases  before  the  expansion  reaches  the  transducers 
is  assumed  to  be  the  projectile  exit  velocity;  with  the  distances  between 
transducers  known  and  rarefaction  wave  time  of  arrival  observed,  the  pro¬ 
pellant  speed  of  sound  can  be  deduced.  With  an  equation-of -state  for  the 
propellant  gas  assumed,  the  propellant  gas  temperature  can  be  estimated. 


For  this  particular  experiment,  a  square— based  plastic  projectile  with 
a  reduced  loading  of  propellant  is  used  in  a  30  mm  WECOM  cannon  with  the 
muzzle  device  attached  as  shown  in  Figure  3.  Five  pressure  transducers  were 
used  in  the  muzzle  device.  The  distances  of  the  gages  from  the  muzzle  are 
given  in  Table  1. 


6.  G.  F.  Widfopf,  J.  C.  Buell,  and  E.  M.  Schmidt,  "Time -Dependent  Near  Field 
Junel1982ake  ?1<M  Sirmlations’  "  AIAA  PaPer  82-00973,  AIAA,  New  York,  NY, 

7.  J.  Huizinga,  G.  Grunwald,  and  N.  French,  "Mobile  ICBM  Transporter  Launcher 
Blast  Loading  and  Hardness  Assessment,  "  Lockheed  Missiles  and  Space 
Company  Report  LMSC-D053108 ,  Huntsville,  AL,  December  1970. 

8.  K.'  S.  Fansler  and  G.  E.  Keller,  "Variation  of  Free-Field  Muzzle  Blast 
with  Propgllant  Type,  "  6th  International  Symposium  on  BaZZ'Lst'icSs 
Orlando,  Fl,  October  1981,  ADPA,  Arlington,  VA. 
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TABLE  1.  GAGE  DISTANCES 


Gage  Number 


Distance  from  Muzzle  (mm) 


1 

2 

3 

4 

5 


7.0 

19.7 
32.3 
45.1 

57.7 


•  H1  f?  propellant  charge  was  designed  to  burn  up  well  before  the  pro- 

J,  I^16  6Xlt®  the  gUn  tUbe  t0  minimize  two-phase  gas  effects.  The  parameters 
of  the  experiment  are  given  below  in  Table  2.  parameters 


TABLE  2.  WEAPON  PARAMETERS 


Bore  diameter,  D 
Bore  and  chamber  volume 
Projectile  mass,  M 

P 

Projectile  velocity,  Vp 
Propellant  mass,  C 
Specific  heat  ratio,  y 
Flame  temperature,  T& 
Specific  force,  F 


30  mm 

u  760  cm3 

112.9  gms 
572  m/s 
12.96  gms 
1.234 
3000  .°K 

1.014*10^  (m/s) 


HI.  NUMERICAL  MODEL 


With  iCieoinVi8^id  m0del  f0r  gUn  blaSt  devel°Ped  fey  Widhopf,  et  al  6  is  used. 
With  this  model,  a  conservative  form  of  the  unsteady  Euler  equations  is  solved 

using  a  one-component  compressible  gas.  These  equations  are: 


(py)t  +  (pvy)  +  (puy)  =  o 
y  x 

(1) 

(Puy)  +  [y(p  +  u2)l  +  (puvy)  =  0 

u  X  y 

/ - N 

CM 

v—/ 

(pvy)t  +  (puvy)x  +  [y (p  +  v2)]y-  p  =  0 

(3) 

(Pey)  +  [puyfe  +  p/p)]  +  [pvy(e  +  p/p)]  =  0 

X  y 

(4) 

P  =  (Y  -  1)  e 

(5) 
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Here  p  is  the  pressure, 
directions  respectively, 
directions  respectively, 
energy. 


x  and  y  are  the  coordinates  in  the  axial  and  radial 
u  and  v  are  the  velocity  components  in  the  x  and  y 
P  is  the  fluid  density  and  e  is  the  specific  internal 


The  simulation  code  utilizes  Godunov's  scheme-9  The  fluxes  at  the  cell 
surfaces  are  determined  by  solving  Reimann's  problem.  Updated  conservation 
values  are  obtained  using  these  obtained  fluxes  in  the  conservation  equations 
with  a  iirst  order  time  differencing.  Surface  properties  are  treated  using 
analytic  expression  for  one-dimensional  wave  reflections  together  with  the 
boundary  condition  that  the  normal  component  of  the  velocity  is  zero. 


his  implementation  of  the  Godunov  method  permits  cells  to  be  sized 
diiferently  to  optimize  resolution  of  flow  structure,  decrease  comortational 
time  and  save  core  memory.  Square  cells  with  dimensions  of  1  mm  were  placed 
throughout  the  bore  from  the  muzzle  to  a  position  somewhat  rearward  of  the 
last  transducer  placed  in  the  gun  barrel  for  the  experiment.  The  cells  in 
the  regions  away  from  the  transducers  are  of  less  interest  and  are  made 
larger.  The  simulation  technique  is  also  modified  to  allow  the  projectile 
base  to  be  initially  at  the  muzzle  plane  position.  This  modification  allows 
a  more  detailed  study  of  the  formation  of  the  expansion  wave. 


he  initial  conditions  outside  the  gun  bore  are  given  as  ambient.  Although 
a  precursor  blast  wave  is  generated  by  the  projectile  pushing  air  out  of  the 
gun  barrel,  this  is  neglected  in  the  present  calculations.  The  in-bore  condi¬ 
tions  are  assumed  to  be  Lagrangian ; 1 0  that  is,  the  velocity  of  the  propellant 
gas  increases  linearly  from  the  value  zero  at  the  breech  to  the  projectile 
velocity  at  the  base  of  the  projectile.  The  propellant  gas  density  is  assumed 

constant  and  the  pressure  increases  from  the  muzzle  to  the  breech  according 
to  the  equation,  b 


P  =  Pm  {(1  -  C/2mp)  [ (L+x) 2/L2] } 


(6) 


where  pm  is  the  pressure  at  the  muzzle  immediately  before  uncorking,  L  is  the 
ength  of  the  barrel  and  x  is  the  axial  coordinate  with  x=0  corresponding  to 
the  muzzle  plane.  The  muzzle  pressure  assumed  for  the  problem  is  the  value 
found  by  a  simple  internal  ballistic  theory.  The  muzzle  pressure  is  83  MPa* 
by  knowing  the  propellant  charge  mass  and  the  bore-chamber  volume,  the  density 
at  uncorking  is  calculated  to  be  170.5  kg/m.3  The  ideal  gas  equation  of  state 
gives  the  pre-uncorking  speed  of  sound  at  the  muzzle,  a,  as  775  m/s. 


9.  S.  K.  Godunov ,  A.  B.  Zabrodin 3  and  G.  P.  Prokopov ,  "A  Computational  Scheme 
for  Two -Dimensional  Nonstationary  Problems  of  Gas  Dynamics  and  Calculation 
of  the  Flow  from  a  Shock  Wave  Approaching  a  Stationary  State U.S.S.E. 
Computational  Mathematics  and  Math  Physics ,  May  1961 ,  pp.  1187-1219. 

10.  J.  Corner ,  Theory  o£  the  Interior  Ballistics  erf  Guns,  John  Wiley ,  NT 

i  nen  ~~  ^  J y  y 
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IV.  RESULTS  AND  COMPARISON  WITH  EXPERIMENT 


Figure  4  shows  a  sequence  of  Mach  contours  generated  by  the  numerical 
simulation  technique.  For  clarity,  Mach  contours  are  displayed  in  a  small 
:  ™nd  2“  SOnJC  line-  The  "1“  line  closest  to  a’doJted  line  is  the 
IZri  f  lu  contours  above  the  projectile  correspond  to  the  recompression 

d^  rlTw6  pr°^llant.8as.Plume*  Initially,  the  wave  front  should  advance 
downward  between  the  projectile  and  the  muzzle  at  a  speed  equal  to  a,  while 
the  wave  front  should  travel  rearward  with  the  speed  a-u.  The  contour  shapes 

^arwLr^r^  T"  lndlcate . that  thia  »  occurring.  As  the  wave  travels 
rearward,  the  contour  lines  straighten  out  rapidly.  Utilizing  the  pre- 

uncor  ing  speed  of  sound  calculated  in  the  previous  section,  the  head  of  the 
wave  should  reach  the  first  transducer  between  the  80th  and  the  90th  steps. 

The  contour  lines  at  the  position  7  mm  from  the  muzzle  are  almost  straight 
a  ese  times.  The  sonic  line  is  also  approaching  the  muzzle  plane  by 
cycles.  Thus,  essentially,  a  flow  pattern  is  established  in  time  to 

obtain  data  for  a  one-dimensional  rarefaction  wave  sweeping  over  all  trans- 
Queers • 


It  is  usually  assumed  that  the  mass  flux  issuing  from  the  muzzle  plane 
has  sonic  conditions  when  performing  a  muzzle  blast  simulation.  This 
assumption  saves  computation  time  and  core  memory.  Examining  the  contour 

in  L8u*e  »  t*le  sonic  line  is  seen  to  be  near  the  muzzle  plane  by 
e  eye  e  of  computation.  At  this  time  , the  muzzle  blast  front  is  still 

quite  close  to  the  muzzle.  Thus,  it  appears  that  the  assumption  is  valid  for 
muzzle  blast  simulations. 

It  is  also  of  interest  to  compare  the  pressure-time  curves  obtained  by 
the  numerical  simulation  technique  with  the  data  obtained  with  the  pressure 
transducers.  Figure  5  shows  the  sequence  of  pressure-time  curves  at  the 
various  transducer  positions  obtained  by  both  experiment  and  calculation. 

n  the  whole,  the  initial  values  for  the  pressure  are  slightly  higher  utilizing 
interior  ballistic  theory  with  full  burn-up  of  the  powder  assumed.  However, 
it  is  likely  that  the  powder  is  not  completely  burnt.  It  is  also  noted  that 
the  pressure  fall-off  is  more  rapid  for  the  experiment  than  for  the  simulated 
results.  The  major  source  of  the  disagreement  is  thought  to  be  caused  by  the 
different  sound  speeds  found  for  experiment  and  interior  ballistics  theory. 

he  experimental  value  was  found  to  be  approximately  840  m/s  whereas  the 
interior  ballistics  value  is  775  m/s  assuming  the  ideal  equation-of-state . 

s sumption  of  the  more  realistic  Abel  equation-of-state  changes  the  calcu¬ 
lated  sound  speed  only  an  insignificant  amount.  Of  course,  faster  sound 
speeds  correspond  to  faster  propellant  gas  emptying  and  therefore  a  faster 
pressure  fall-off.  Studies  indicate2’3  that  the  propellant  gas  near  the 
base  of  the  projectile  could  have  a  higher  temperature  than  in  the  region 
toward  the  breech.  The  lack  of  better  agreement  may  also  be  caused  partially 
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nLnhe  1?rg?  artifi?ial  viscosity  implicit  in  the  first  order  Godunov  tech- 
q  e.  mplementation  of  second  order  Godunov  schemes11* 12  might  produce 
a  crisper  more  definitive  waveform  and  better  agreement  with  experiment. 

f.  V-  6  one_dimen8ional  assumption  is  valid  for  the  projectile  in  the 

tie  result  wienM  there  ouSht  to  be  good  agreement  between 

tne  results  when  the  projectile  was  included  and  when  the  projectile  was  not 

tUUvte  dt  “  th?  «»  fl°»  field  within  the  Ire  it  eesen- 

tially  one-dimensional  immediately  after  firing  and  remaina  ao.  Figure  6 

^  the^rSld0"  b?JWee"  the  Pr“8Ure  °btai”e'i  ”Uh  “d  »  Projectile 

llnAiti  J  Ui  lhXS  se(luence  was  run  using  slightly  lower  muzzle  exit 
calcnl^ t  a  ^an,  the  former  sequence.  As  might  be  expected,  the  pressures 
caicuiated  at  the  transducer  are  at  first  identical.  Later,  the  effect  of 

“ST?  at  and  near  the  muzzle  is  seen  to  cause  the  pressures  to  be 
slightly  higher.  Nevertheless,  the  general  shape  of  the  curve  is  unaltered 

identical 8PeedS  dedUCed  fr°m  the  two  cases  would  be  expected  to  be  nearly 


VI.  CONCLUSIONS 

For  all  practical  acoustic  thermometry  experiments  of  interest,  the 
one  dimensional  assumption  for  the  rarefaction  wave  appears  valid.  Compari- 
be tween  simulation  and  experiment  yields  fair  agreement.  The  maior 
source  o  isagreement  stems  from  the  higher  sound  speeds  obtained  by  experi¬ 
ment.  The  reason  for  the  higher  sound  speeds  is  not  known  although  there 

mod!Ilbeen  S°m!  aXp^anatlons  for  this  behavior.3  Of  course,  if  the  numerical 
odel  were  modified  to  take  into  account  the  possibility  that  higher  tempera- 

much  °btained  near.the  muzzle,  the  results  from  the  numerical  model  would 

more  closely  approximate  the  experimental  results.  Most  previous  experi- 
menters  have  obtained  lower  sound  speeds, but  some  of  the  results  could  be 
attributed  to  incomplete  burning. 


11  ‘  Bram  van  Leer ,  "Towards  the  Ultimate  Conservation  Difference  Scheme.  V. 

A  Second  Order  Sequel  to  Godunov's  Method,"  J.  Comp.  Thus.  (32) 
pp.  101-136 J  1979. 

12.  S.  Eidelman ,  P.  Colella ,  and  R.  P.  Shreeve,  "Application  of  the  Godunov 
Method  and  Higher  Order  Extensions  of  the  Godunov  Method  to  Cascade 
Flow  Modeling 3  "  AIAA-83-1941-cPl  AIAA  6th  Computational  Fluid  Dynamics 
Conference ,  Danvers ,  MA,  July  13-15 3  1983. 
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MUZZLE  EXPANSION 


Figure  1.  In-Bore  Wave  Diagram 
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Figure  2.  Pressure-Time  Trace  at  Transducer  Station 
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Figure  3.  Experimental  Setup 


(a) 


(b) 


Figure  4.  Mach  Contours  -  In-Bore  and  Around  Muzzle  (continued) 
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(a) 


19.7  mm  Y=  16.0  mm 


(b) 


Figure  5.  Pressure-Time  Curves  -  Simulation  vs  Experiment  (continued) 
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Figure  5.  Pressure-Time  Curves  -  Simulation  vs  Experiment  (continued) 
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x=  -57.7  mm 


Y  —  16.0  mm 


Iigure  5.  Pressure-Time  Curves  -  Simulation  vs  Experiment  (continued) 
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